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A one-dimensional quantum N-body system of either fermions or bosons with 
SU(n) "spins" (or colors in particle physics language) interacting via inverse-square 
exchange is presented in this article. A class of eigenstates of both the contin- 
uum and lattice version of the model Hamiltonians is constructed in terms of the 
Jastrow-product type wave function. The class of states we construct in this paper 
corresponds to the ground state and the low energy excitations of the model that can 
be described by the effective harmonic fluid Hamiltonian. By expanding the energy 
about the ground state we find the harmonic fluid parameters (i.e. the charge, spin 
velocities, etc.), explicitly. The correlation exponent and the compressibility are also 
found. As expected the general harmonic relation(i.e. v$ = (vjsrvj) 1 ^ 2 ) is satisfied 
among the charge and spin velocities. 
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I. INTRODUCTION 



In this paper we report an exact treatment of a one-dimensional system of particles 
with SU(n) "spins" interacting with the inverse- square, two-body exchange in continuum 
space. In particular, we show that a class of eigenstates of the Hamiltonian is given by the 
Jastrow-product type of wave function. We also show that this model can be put on the 
lattice for all the positive integer values of the dimensionless interaction parameter which is 



defined in Eq. ( |2.1|) . The continuum version is a direct generalization of Sutherland's model 
of interacting spinless fermions (or hard core bosons) |l| . The lattice version corresponds to 
the SU(n) t-J model which, at n — 1, maps to the Haldane-Shastry model p. 

Fig. 1 shows how various models with the inverse square exchange have evolved from the 
original Sutherland's Model. In the figure the models that already exist are in rectangular 
boxes and the new models that we present in this paper are in ovals. The solid arrows 
mean "generalized from" or "evolved from". Two equivalent models that are represented 
in different ways are connected by a dashed arrow. We emphasize that while the bosonic 
t-J model (A odd) is a direct generalization of the bosonic supersymmetric case (A = 1), 
the fermionic t-J model (A even) is not directly related to the fermionic supersymmetric t-J 
model that Kuramoto and Yokoyama H found. 

II. MULTICOMPONENT GENERALIZATION OF SUTHERLAND'S MODEL 

For convenience we consider a periodic system, where the Hamiltonian in units of h 2 /m 
is given by 

where A is the dimensionless interaction parameter. P-j is an operator that exchanges particle 
spins at Xi and Xj, and d(x) = (L/ir) \ sm(irx/L)\. d(xi — Xj) is the chord distance between 
particles at Xi and Xj on a circle with circumference L. If all particles have the same spin, 
this model reduces to the system of spinless particles studied by Sutherland. Note that 



Sutherland's coupling parameter A' corresponds to A — 1 in our notation. For example, the 
spinless free Fermi gas corresponds to A' = 1, but in our notation A = 0. 

It is easily found that the wave function \1/ must vanish j\ ^~ (1^2 I ) 

Xi — Xj — > in the case of symmetric (antisymmetric) spin configuration of the two particles 
at X{ and Xj. And, if < A < 1 the effective interaction strength is attractive for the 
antisymmetric spin configuration and there is some ambiguity in this it must be 

further specified whether \1/ vanishes particles approach. If 

we choose the first boundary condition the free fermion limit is obtained as A — > 0. 

There are two possible interpretations for this system. The first is to consider the system 
as a ring embedded in a plane. Hence, even though the particles are constrained to move 
in one-dimension, the interaction is two-dimensional in nature. The other is to regard 
the system as strictly one- dimensional by taking l/d(xi the effective interaction 

after summing the pairwise interaction of the particles around the ring infinite times. The 
following identity, 

+oo ^ ^ 

^(x + nL)^^ (2 ' 2) 

has first been used by Sutherland [0 to show the validity of the latter interpretation. 

In analogy with the states previously constructed for the SU(2) spin chain in ref. 0, we 
propose the following Jastrow-product type wave functions for our Hamiltonian: 

*({M) = *on4 CTfc , (2-3) 

k 

where 

^0 = II 0nm 

n>m 

<Pnm = \zn - z m \ x ~ x (z n - z m ) x+5,T " am exp(i^7Tsgn(a n - a m )). (2.4) 

Here, 5 is the Kronecker delta function and z n = exp(2nix n /L). a n is the ordered spin index 
and J rj the global current of particles with a spin. We take J a to be an integer and will 
discuss the restrictions on the allowed values of J n later in this article. Note also that the 



wave function with A = x = is the Slater determinant that corresponds to the states of 
free SU(n) fermions. 

The symmetry of the wave function with respect to the exchange of particles is given as 

*(..., Zi<Ti, ZjCTj, ...) = (-l) x+1 tf (. . . , Zj<Tj, Zi<Ti, ...). (2.5) 

Hence, for bosons(fermions), x — l(x — 0). 

We write the total Hamiltonian as H = H° + H 1 + H 2 where H°, H 1 and H 2 are the 
kinetic, potential and spin exchange Hamiltonian, respectively. We will show that each 
operator acting on the wave function gives two types of terms, "wanted" and "unwanted". 
"Wanted" terms are defined to be the terms that depend only on the global variables such 
as the total number of particles, J a , etc. "Unwanted" terms explicitly depend on the local 
variables such as Zj. Since the eigenenergy should depend only on the global variables, the 
"unwanted" terms for if , if 1 and H 2 should cancel or combine to give "wanted" terms. 

We first examine if acting on the wave function. We define the following derivatives: 

^^ = -^-^--%L_ (2.6) 

€ij = d Zj (pij (2.7) 



*? = z rip- = z i E (2-8) 



i&3) 

.f-^ = ^E^ + (^) 2 . (2.9) 

In terms of these derivatives we find that H ^ = 2(n/L) 2 Ejfaf } + (2J ffj + l)^ + 
After some algebra one can show that if acting on the wave function gives "wanted" (W K ) 
and "unwanted" (U K ) terms. We write these terms as follow: 

L, 



H°* = 2l?r) (W K + U K )V, (2.10) 



where 



W K = ^N(N 2 - 1) + i(A + 1) ]T M -( M - " !) + I E M -( M - " !)( M - " 2 ) 



- ]T M a J a (XN M a J a A + l) + E{x), (2.11) 

where £(x) = x 2 iV(jV - l) 2 /4 + (x/2)(iV - 1) £ CT M a (M a - 1) + x(JV - 1) E CT M CT J CT . M a is 
the total number of a particles in the system. We can now anticipate Ur to be canceled or 
combined with the "unwanted" terms from H 1 and H 2 to give local variable independent 
terms. 

H 1 acts trivially on \I/ to give A 2 YU<j ~^/d{xi — Xj) 2 . The action of H 2 on the wave 
function is less trivial. We need to evaluate the following expression: 

p w = i E ir^z = -^ttE j-^nW}), (2.i3) 



Zi Zi 



where {za 1 } is a configuration with cxj and <ij exchanged with respect to {za}. Similarly, 
we let {z'a} be equal to {za} with Z{ and Zj exchanged. Using the identity, \P({zcr}) = 
(— l) x+1 ^({z'cr'}), we may then rewrite the exchange operation as 

z iZj ^>{{za'}) 



{M} 



-ir +1 E 



^(Zi-zjy^z'a'}) 



= e ^^(-i) 5 -^ - n ZJL -^ ) (2-i4) 

i^y l z « z j) \ z j ) k^ij \ Zk z i ) 

The action of P a is same for both boson and fermion cases and is independent of the 
interaction parameter A. 

We use the following identity, 

s m /, ,v f/wy s(n) + /v-£\ 5 S( _„)\ , (2 . 15 ) 



9=1 

to rewrite Eq. (|2.14[) as 



Q 



^/i = ^o + E E p f- ( 2 - 16 ) 

<j,<j' q=l 



9 in Eq. (|2.15| ) is the step function with 8(0) = 1/2, and P an d P™' are given by 
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/ \ ^ <7 j (7 k ^ a i a k 



pf = 2 e Tr^h^ wi - <w n 



ijkj \ Z i ~ Z j) k^ij \ Z k ~ Z jJ 

x(^)e(W„)(^) a . (2.18) 

It is useful to separately consider the terms with q — 0, l and the terms with q > 2. For 
the terms with q = 0, l we introduce two sets of site indices {a} and {f3}. The set {a} ({/?}) 
includes all the locations of particles with the spin cr(cr'). Using the identity, 

(^1^'" = l - ^ *Z?L + 5 aak ^i- for a ^ a', (2.19) 

the products in P{m} m ay be expanded and the typical terms in P{m} can be simplified 
using the following two theorems. 

Theorem 1: Let {a} and {j3} be two sets of distinct integers between 1 and N, and let 
< q < 1 and A = (1 - 5 a y)5 aai 5 a(7a - ■ ■ 6 <T<Tan 6 <T > tT .6 IT > <T £ ■ ■ <5 CT ' CT/3m - Then, 

M CT -1 M a ,-1 (l\m ? >?/ r _ _ \n+m-2+<j 

y^ y^ y^ y^ I L ) flfj ^j/ ^ 

- Ef=T (MCT ' MCT,) (M CT - jfe)(M^ - fc) for g = , 

_ E mn(M M{Ma _ k) forg=1 . 

Theorem 2: For g > 2 the following identity will hold for M a > M a >, 

e - f^^V V- www, n f^^V CT ' CTfc 

i^j % V Z J I k^ij \ Z k Z j / 

M a , for q = 2, 

for 2 < q < M a - M c , + 1. 

Proofs of these theorems are given in Appendix A and B. 

After reorganizing some terms we obtain the following results: 

P {M} = W P + Up, (2.20) 

where 



Wp 



--Y l M a (M a -l)(M a -2)--Y t M a ,(M a ,-l)(3M a -M a ,-l), 



(2.21) 



Up 




(2.22) 



Because of symmetry we can choose M\ > M 2 > . . . > M n and < J a > — J a < M a — M a > + 1 
for a' > o without loss of generality. And, one can easily check that Uk, U p and the potential 
energy term combine to give local-variable independent terms. 

Before we give the full expressions for the eigenenergy, we discuss the allowed values of 
the integer currents. Theorem 2 gives a simple rule for selecting allowed currents. If we 
choose -1 < J a < — J a < M a — M a > + 1 for M a > M a ,, then P q = for q > 2. Otherwise, the 
energy will not in general be independent of the local variables and the corresponding wave 
function will not be an eigenfunction. Pictorial illustration of the allowed currents is shown 
in Fig. |3|. Each row of M v boxes represents M v "particles" of same spin. A single box gives a 
unit of current. Fig. |](a) is the reference state with no current. If the current is j > 0(< 0) 
for a row, the row is moved j boxes to the right (left) as shown in Fig. 0(b). In order for a 
state to be an eigenstate, the following condition must hold true: For a given pair of rows 
of boxes all except the first and the last boxes in the row with smaller number of boxes must 
be positioned within the other row. Fig. |3|(c), for example, cannot be an eigenstate because 
the last two boxes in the second row are not within the first row. 

Finally, the eigenenergies of the model can be given as follow: 



where E\ and E 2 are energies due to one- and two-spin interactions, respectively, and are 
given by 




(2.23) 



E 



~ 12 



N(N 2 - 1) + M a (M a - 1) + Ul - \)Y,M a (M a - l)(2M a - 1) 
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+ £j ff M CT (M ff + J*-l) (2.24) 
E 2 = ~h E M °\ M °> - 1)(3M CT - M a , - 1) - A E M CT ,(M CT - M^J^ - J CT ) 

d <T<er' a<a' 

+ M a ,{J a/ -J a ) 2 (2.25) 

<T<0-' 

The ground state for fermions (bosons) is given by the following two conditions: (i)M a > = 
M a = M and {ii)J a , = J a = J = — (M - l)/2, -(N + M - l)/2 for fermions and bosons, 
respectively. The currents should be integers. If the condition (ii) does not give an integer 
current the ground state is degenerate and the ground state currents are two integers closest 
to the the half odd integer value. 



III. SU(N) T-J MODEL 

In ref. 0, Sutherland's model for a system of spinless particles in continuum space was 
extended to the lattice. Using the similar procedures we put our model on the lattice. We 
propose the following Hamiltonian: 

I is a positive integer. N a is the total number of sites and N the total number of particles. 
d(n) = (N a /ir) \ sin(n)|. V- hops a particle at jth site to (j +n)th (mod N a ) site, rij is equal 
to one(zero) if jth site is occupied (empty). is, as before, the spin exchange operator. V 
is the projection operator that insures the absence of multiply occupied sites. 

In analogy with ref. we propose the following eigenstate for the model Hamiltonian, 

IV})=E«IW>- (3-2) 
W 

Since the total number of each spin is a good quantum number we classify the eigen- 
states into sectors labeled by {M} = (Mi, M2, ■ ■ ■ , M n ), where M v is the number of 
z/th spin and YZ=\ M u = N. We also represent the particle configuration as {za} = 
(zicri, . . . , ZiOi, . . . , Zjdj, . . . , z N a N ), where and <jj are the location and spin of the zth 
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particle. The sum in Eq. ( [3.2|) is over NJ/ \Mh\M\\ ■ ■ ■ M n \) distinct spin configurations of 
the sector {M} where N a and are the total number of sites and the number of holes, 
respectively. The function </> in Eq. ( |3.2| ) is given by the following Jastrow-product, 

0({4IM) = II(* -z^^e^^-^Hz^. (3.3) 

i<j k 

Considering the symmetry of the wave function we take odd (even) I to be bosonic (fermionic) 
state. 

As in the continuum case we break up the total Hamiltonian into the kinetic, potential, 
and spin exchange parts (H®, H\ and iff, respectively). The actions of H\ and H\ are 
same as the corresponding operators in the continuum case. For H\ we have the following 
relations: 



<H U > 



{ziai, . . . , z N a N \ 


H°\ 




(zioi, . . . , z N a N \ 




{M}) 



=4E e ^(i-^a-o- 1 n ^— ^ (3-4) 

i=l n=l jtyi) \ Zl Z i / 

where \z\Gi, . . . , zn&n) is one of the basis states of the sector |\I/{m})- In order to evaluate 
< H\ > we need the following theorem (a slightly different version of this theorem can be 
found in ref. 0). 

Theorem 3: Let J and p be non-negative integers and z = exp(2ni / N a ) . Then, 
EV J (1 - z- n )-\l - z n f- x = - J{Na - J) for p = 0, 

n=l 12 2 

= —J H — torp=l, 

= +1 for p = 2, < J < N a - 2 
= -(JV a -l) for p = 2, J = N a — 1 

m is a positive integer, and the restriction on the current is < J < iV a — 1 unless specified 
otherwise. 
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If we take I > 0, we can multiply out the product in Eq. (|3.4j) . The resulting expression 
will be a sum of terms containing a factor (1 — z n ) p . Since the maximum value of p is 
l(N — 1) + M CT — 1, if we choose < J a < N a — l(N — 1) — M a + 1, we may ignore terms 
with p > 3. We combine the remaining terms with those from H\ and H\ and obtain the 
following eigenenergies for the lattice model: 

E L = E{ + £ 2 L (3.5) 

where 

El = ±N(N 2 a - 1) + ^N(N - 1)(JV - 2) - l -(N a - l)N(N - 1) 

- {^y^ ~ l ( N ~ x )) E M -( M - - 1) + |(i - E M -( M - - - 2 ) 

- £ M CT J CT (iV a - ZJV - M CT - J CT + Z + 1), (3.6) 

(J 

£ 2 L = E 2 , (3.7) 



where -E2 is given by Eq. ( 2.25 ). We emphasize at this point that the above expressions for 
the eigenenergies for the lattice model are valid only for IN + M a < N a + I + 1. 

For / = 1 (bosonic supersymmetric case) we can find another expression for the energy 
for iV + M a > N a + 2. We proceed by rewriting the kinetic energy term as 

where o~h represents the hole. The above expression is almost identical with the expression 
for the spin exchange operator and can easily be evaluated using Theorem 1 and 2. The 
energy then is given by 

E(l = l) = ^(E k (l = l) + W P ), (3.9) 



where 



Ek (l = 1) = _ Mfe ^ jj _ \ Mh (M h - l){M h -2) + M h Y J J*(J* ~ M a - 1) 
6 3 V 

+ (E M -( M - - X ) + M ^ - !)) - - 1)[3JV - n(M ft + 1)] 

+ M h N + Ea " 1)[2AW4 - M h (Af fc + 1)], (3.10) 
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where the currents are restricted to —{M a — Mh) < J a < 0. Wp in Eq. ( |3.9| ) is given by Eq. 



For / = the wave function in Eq. ( |3.2| ) is that of Gutzwiller projected free fermions 
and no longer an eigenstate of the Hamiltonian in Eq. ( |3.1| ). However, the wave function is 
an eigenstate of the following Hamiltonian 

ff = P (_ E _5_ +E (3zWU ,3.11) 

This Hamiltonian is SU(n) generalization of the SU(2) fermionic supersymmetric t-J model 
which Kuramoto and Yokoyama || solved. If we employ the similar technique used for I = 1 
case we may find the eigenenergy given by Eq. ( |3.9|) with Ek(l = 1) replaced with Ek(l = 0) 
which is given by 

E k (l = 0) = _^lZ22 + h Na _ i) £ Ma (M a - 1) + M h Ja{M a + J a - 1). (3.12) 

£7 

The currents are integers restricted to — M CT < J CT < 0. The ground state energy is obtained 
if J rj = — (M a — l)/2 and M a odd. If M a are not all odd, the ground state is degenerate. 

IV. SU(N) SPIN CHAIN 

If we take the spin exchange Hamiltonian as our full Hamiltonian for a system of spins 
on a lattice ring, we obtain the SU(n) spin chain model. This is equivalent to taking the 
"half filling" limit (N —>■ N a ) of the lattice t-J model. In this model the eigenenergies for the 
boson and fermion systems are identical. Since every site is occupied with a spin, Up in Eq. 



(p|) can now be summed and is equal to -A^ a (iV a 2 - 1) /12 + {(N a - 1) /2} £ CT M a (M a - 1 



Hence, Up is now local variable independent, and the eigenenergy is given by Eq. ( |2.20|) . 



V. THE HARMONIC FLUID DESCRIPTION 



We now turn our attention to the low energy excitation spectrum of our model in the 
thermodynamic limit. The low energy excitation spectrum of a one- dimensional quantum 

11 



fluid may be described by a sum of two independent harmonic fluid Hamiltonians, one for 
the charge and the other for the spin. This is a slight generalization of the spinless particle 
system in ref. 0. In general the effective Hamiltonian can be written as follow: 

H = J2jdx (A aa ,Il a (x)n a ,(x) + B^V^^V^ix)) (5.1) 

Htr(x) is the local density fluctuation field as in ref. ||. 4> a (x) is the canonical conjugate 
field to Il(x), [(p a (x), ILy(x')] = i5 aa >8(x — x'). Because of the SU[n) symmetry we have 
A aa i — ai + ao5 aa i and B aa > = b\ + b 5 aa ' . We now express the Hamiltonian in terms of the 
Fourier transformed fields, U ak and <fi a }.. 

^ = EE A aa'^ak^a'-k + B a y k 2 <f) uk <f) af _ k . (5.2) 
k <t,ct' 

We now construct normal mode fields as U k = J2a a a^k and (fi k = J^aK^^k- If we choose 
J2aK a a = <Wj then [0^,ll^' fc ,] = iduk'Svu'- We have the following equations of motion: 
[H, [H,Ul\] = -(v u k) 2 W k , [H, [H,<f>l)) = -(v u kf<p u k . For H-field we obtain the following 
equation, 

(a b - (v l/ ) 2 )a^ + (ai&o + a &i + nai&i) = (5.3) 

The same equation with a v a replaced with b v a holds for 0-field. There are only two possible 
values for (v u ) 2 : a b and (a + nai)(b + nbi). The first value corresponds to the case 
J^pttp — and the second J2/3 a p 0- Hence, the first would be the spin velocity and the 
second the charge velocity. The Hamiltonian can now be written as 



k 



n— i 

(oq + n ai )U c k WL k + (b + nb x )k 2 <\> c k <\> c _ k + £ {a Q n£n s _ fc + b k 2 <f> s k (f> s _ k } 



n-l 

[ (5.4) 

s=l 



We have one charge mode and n—1 spin modes; and because of the SU(n) symmetry all the 
spin modes have same velocity. We define the following velocities associated with n- and 
0-fields for the charge(denoted as c) and spin(denoted as s) modes: v% = a + na 1; v 1 } = 
bo + nbi, v s N = a , Vj = b , . The charge and spin velocities then are v c = (v c N v c j) l l 2 and 

V s = {v S N V S j) 1/2 . 
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Since the eigenstates for our model corresponds to the low energy excitations of the 
model, we can easily obtain the harmonic fluid parameters by examining the energy expanded 
about the ground state. By Taylor expansion we obtain the energy for the continuum model 
as follow, 

E = eoN + fiN+ ^l^A J2 l -A ary , AM a AM a , + 2B aa , AJ a AJ a , (5.5) 

771 Lj , Z 

where 

eo .(^A + lV, (5.6) 



6m V 7i/ 
(nhpo) 2 ( l\ 2 

=a(a + -) + (a+£)<W, (5.8) 
B aa , = --+(\ + -)5 aa ,, (5.9) 



where p — N/L. We can now just read off the velocities, 



m V nj 


(5.10) 


v c N = n A + - 

m V nj 


(5.11) 


c 1Ttip 1 
Vj = 

m n 


(5.12) 


V c = V s 


(5.13) 



The charge and spin velocities are same for all A and n. As expected from the singlet nature of 
the ground state the ratio Vj / v s N does not get renormalized. The coefficient Vj is independent 
of the interaction term, as a consequence of Galilean invariance, but v c N gets renormalized due 
to non-linear interaction terms not included in the linear form of the harmonic Hamiltonian. 
The renormalization coupling constant for the charge is exp {-2<f) -> {v c N /v c J ) 1 ' 2 = n\ + 1 
in the limit of long wavelength (i.e. k — > 0). Because of the scale invariance of the model 
the dimensionless coupling constant if is independent of the particle density. 

The compressibility per particle is (n 2 U 2 pl/m)(\ + 1/n) 2 . We also find the chemical 
potential and the ground state energy to be m*v 2 F /2 and (N/3)(m*Vp/2), respectively. The 
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chemical potential (or the Fermi energy at zero temperature) and the ground state energy of 
this one-dimensional system is that of free fermions with the renormalized mass per particle 
m* = m/(n\ + l) 2 . The Fermi velocity is given by Vp = nhpo/(nm*). 

The low energy excitation of the lattice version is given by Eq. (|5.5|) as well if we set 
m = h = 1 and replace L, e and p with N a , — 7r 2 /3 + 2e , and — 7r 2 /3 + 2p, respectively. The 
crucial difference is that this expression for the lattice model is true only for po(= N/N a ) < 
pmax^_ (7_|_i Hence, the charge velocity is linear in p only up to p™ ax . We expect v c 
to vanish as po — > 1 due to the metal-insulator transition at the density. v c , therefore, should 
exhibit non-analyticity at p™ ax . This behavior is attributed to the long-range interaction in 
our model. In ref. JiJ this type of non-analyticity in the spinon velocity was observed for 
the SU(2) spin chain model. In the nearest neighbor interaction models the sharp change 
in the charge velocity is smoothed out. 

For I = 1 the energy for po > (1 + l/n)^ 1 is explicitly found and the harmonic fluid 



Hamiltonian is given by Eq. ( |5.5| ) with A aa > = (1 — p )(n + 2) — p /n + 5 aa i and B a(T > = 
—po/n + 5 ua i. Therefore, v c = n(l — p )(n — 1) and v s = ir. The spin velocity is independent 
of the density in this region. Similarly, for / = we have A aa i = B aa i = —po/n + <5 CTCT / for 
< Po < 1- Thus, v c = 7r(l — po) and v s = tt. This agrees with the results obtained by 
Kuramoto and Yokoyama || . However, it is unexpected that the charge and spin velocities 
are independent of the number of spin species. Fig. 4 shows the charge and spin velocities 
for the lattice models with various interaction parameters. 



VI. CONCLUSION 

We have shown in this paper that the Sutherland's model of spinless particles inter- 
acting with inverse square exchange can be generalized to the multicomponent system of 
particles. We further generalized the continuum model to lattice for the integer values of 
the dimensionless interaction parameter. We explicitly constructed a class of eigenstates 
for these models and calculated the corresponding eigenenergies. We claimed in this article 
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that the class of eigenstates correspond to the ground state and the low energy excitations 
of the model. We checked this by the exact numerical diagonalization of the lattice version 
of Hamiltonian for small systems. The systematic construction of the other eigenstates, 
thermodynamics, and the general n-point correlation functions for these models are yet to 
be found. As J. A. Minahan showed in his preprint J7| the Hamiltonian of the type given by 
Eq. (|2.1| ) also arises in a matrix model as a representation of one dimensional open string 
theory. 

This work is supported in part by NSF grant No. DMR-91-96212. 

Note added: After the completion of this paper we received a preprint from N. Kawakami 



who reports similar results. And, the Hamiltonian given by Eq. ( [2.1| ) is also independently 
identified by A. P. Polychronakos ||. 



APPENDIX A: PROOF OF THEOREM 1 

We first discuss a diagrammatic way of writing the exchange operation. A labeled di- 
agram is shown in Fig. |2|. The amplitude of the labeled diagram shown in Fig. |2| can be 
evaluated using the following rules: i) the dashed line connecting the indices i and j gets a 
factor ZiZj(zi — Zj)~ 2 { \ — 5 a ^i)5 uai 5 a(T:j1 ii) each solid line connecting the indices i and gets 
a factor (zi — Zj)(z ak — z i )~ 1 5 a(Tak , iii) each solid line connecting the indices j and (3k gets a 
factor —{zi — Zj)(z/3 k — Zj)~ 1 5 aapk . Weight of the dashed line vanishes whenever the spins at 
the sites i and j are same. On the other hand, weights of the solid lines that are connected 
to vanish whenever the spins at the sites {a} ({/?}) and i(J) are different. Hence, we 
may think of the diagram as two-spin interaction diagram. Since the Hamiltonian has only 
the two-spin exchange operator Py, three-spin and higher number of spin interactions are 
missing. 

The amplitude of the labeled diagram is undesirable since it depends on the local vari- 
ables. A more desirable, local variable independent diagram (unlabeled diagram) can be 
obtained by summing {n + l)(m + 1) cyclically permuted labeled diagrams. This fact will 
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be proved in this appendix. 

We note that the indices, {i, cui, . . . , a n } and {j, /?!,..., j3 m }, are dummy indices. Conse- 
quently, the sum will be invariant under any permutation of the indices. We further observe 
that A is invariant under any permutation of the indices except when there is interchange 
between the two sets, {i, a 1; . . . , a n } and {j, Pi, ... , /3 m }. On the other hand, the pre-factor 
to A remains invariant with respect to the permutations of indices {a} and {/?}. It is conve- 
nient to consider only the cyclic permutations which give (n + l)(m+ 1) distinct pre-factors. 
We, then, sum the (n + l)(m + 1) factors as follow: 



pcrcr 

nm 



E E 

¥j Mm 



Z;Z 



1-9, 



\n-\-m— 2-\-q 



! ml ( Zi - z ai ) ■ ■ ■ (zi - z an )(zj - Zfi) ■ ■ ■ (zj - z p J 



n 



■iy 



X 



i 

A, 



(Al) 



where {7"} = {i, a x , . . . , a n } and {7^} = {j,p 1 ,..., f3 m }. 



We let D r 



Yl^i'iz — z 0). Using the binomial theorem we 



expand the numerator in Eq. (|A1|) . We then multiply and divide the expression by D a Dp 
and obtain the following. 

_^m-l+ q -s / n + m _ 2 + q \V°Vp r ' 



n+m— 1+q 

= E E E ■ 

i+j {a},{/3} 8=1 



rn+rn—s 



n + 1)! (m + 1)! 



1 



D a Df 



-A, 



(A2) 



where 



V s 



Zi z. 



ra-1 



v n-l 



Oil 



... 2: 



(A3) 



and Vg is defined the same way. We further note that V^ = 0if0<s<n — 1 and Va = 
if < /c < m — 1. Therefore, if < q < 1, the only non-zero contribution in ( |A2|) is the 



term with s — n. Since VT = -D„ and V, 



P 



Dp, we obtain 
16 



nm 



(n + l)!(m + l)!V n-1 J^ { J^ } 

'n + m - 2 + g\ / M a \ ( M a , N 
. n — 1 )\n + l)\m+l. 



:-i) 



n+m—l+q 



(A4) 



Evaluating S{a},{/3} A in (|A4|) is same as calculating the total number of ways of putting 
n + 1 out of M a blue balls in a box and m + 1 out of M a > red balls in another box. The sum, 
therefore, is equal to M^M^ - 1) • • • (M a - n)M a ,{M a , -!)■■■ {M a > - m). 
It is straightforward to prove the following useful identities: 

'n — l\ (m — 1 + q\ 



n + m — 2 + q 
k n-1 , 



E 

s=0 

s 

E( 

fc=0 

£(-l) fc ( S + l-£;) 

fc=0 



m + q 
k i 



m + 1 + q 



(A5) 



(A6) 



Mcr-l 

E (-!) r 

m=l— 5 



v m + 1, 



m + 1 



-M CT + 1 if k = 
-M„ if = 1 



> for q = 0, 



for 2 < A; < M CT 
1 for fc = 

for q — 1, 

for 1 < A; < M CT " 



(A7) 



Using (|A5| ) we may sum the terms depending on n and m in separately. Using (A6) 
and ( |A7| ) the following relation can easily be obtained: 



pc 



E M 4 „(M CT ,A V ) (M(j _ k - m ^ _ k) ioiq = Q 
Y%li {Ma,M °'\M„-k) forg=l. 



(A8) 
Q.E.D. 



APPENDIX B: PROOF OF THEOREM 2 



We replace the product over k with two new sets of integer indices, {a} and {j3} and 
rewrite the expression in Theorem 2 as, 
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g ' £j | Zi Zj ^ £ ai Z a2 Zj Za^p—i Zj 



^(^-^(^-l)!^ V h 



E (* - *i) r fr{{zp}) ) A. (Bl) 

r=0 



Here, A is same as in Theorem 1 . The factor in the curly bracket in (|B1|) is obtained by 
writing the terms like (zp l — Zi) / (zp t — Zj) as {1 — (zi — Zj)/ {zp l — Zj)} and by multiplying out 
all the factors that depend on {/?}. f r in ( |B1|) is some function of {zp}, however, f = 1. 

As in Appendix A we rewrite the expression by summing over the factors obtained by 
the cyclic permutations of the indices {i, a%, . . . , qm„-i} and obtain 

Q = 5 S w-d^ - - *> ■ ■ • - < B2 > 

where W^f" 

is given by the Vandemonde determinant whose last row is modified to z^ t (2^ — 
z j) q+r ~ 3 j ' z q f l i = 0> a ii • • • > a Ma- 1}- ^ * s straightforward to show that 

f 7 Vl bigffil ifg = 2andr = 

^yMc _ I l^i z j)( z a 1 Zj ) — (.Zat Ma -1 ZjJ (B3) 

[ ifO<g + r- 3<M CT -3 
E E A = MJM CT ,! (B4) 

#j {a},{/3} 

The only non-zero contribution to Q is given by the condition q = 2 and r = 0. Since 
< r < M CT / - 1, the sufficient condition for W™ a = is 3 < q < M a - M a , + 1. Therefore, 
we have for M CT > M CT / 

M CT / for q = 2, 

Q={ (B5) 
for 3 < q < M a - M a , + 1. 

Q.E.D. 
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FIGURES 

FIG. 1. Evolution of the inverse square exchange models both in the continuum space and 
lattice. Models that already exist are in rectangular boxes and new models that are presented in 
this paper are in ovals. The models enclosed in the dotted box are continuum version and the rest 
of the models are lattice version. The solid arrows mean "evolved from" or "generalized from". 
The dashed arrows connect two equivalent models in disguise. 

FIG. 2. Labeled diagram for n u-spins and m cr'-spins. 

FIG. 3. Illustration of allowed values of the currents for the sector {14,12,10,6,4,2}. (a) 
Allowed. All currents are zero, (b) Allowed. Ground state for the sector, (c) Not allowed. 

FIG. 4. Charge and spin velocities of the lattice models, (a) Fermionic supersymmetric t-J 
model. The velocities are independent of the number of spin species, (b) Bosonic supersymmetric 
t-J model. Up to the density n/(l + n), v c and v s are equal, (c) Bosonic and fermionic t-J models 
with I > 2. The charge and spin velocities are same up to the density (I + l/n) _1 . Beyond the 
density the ground state is no longer the Jastrow-product type. 
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Fig. 3 
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